Wave function collapses in a single spin magnetic resonance force microscopy 



O 

o 

Oh- 
< 

o 

(N 



> 

in 
o 

o 

o 

^: 
Oh: 

^— > ' 
G ■ 

3 : 
cr. 



G.P. Berman, 1 F. Borgonovi, 2,3 and V.I. Tsifrinovich 4 

1 Theoretical Division, MS B213, Los Alamos National Laboratory, Los Alamos, NM 87545 
Dipartimento di Matematica e Fisica, Universita Cattolica, via Musei 41, 25121 Brescia, Italy 

s INFM, Unita di Brescia and INFN, Sezione di Pavia, Italy 
4 IDS Department, Polytechnic University, Six Metrotech Center, Brooklyn, New York 11201 

(Dated: February 1, 2008) 

We study the effects of wave function collapses in the oscillating cantilever driven adiabatic rever- 
sals (OSCAR) magnetic resonance force microscopy (MRFM) technique. The quantum dynamics 
of the cantilever tip (CT) and the spin is analyzed and simulated taking into account the magnetic 
noise on the spin. The deviation of the spin from the direction of the effective magnetic field causes 
a measurable shift of the frequency of the CT oscillations. We show that the experimental study 
of this shift can reveal the information about the average time interval between the consecutive 
collapses of the wave function. 

PACS numbers: 



The oscillating cantilever driven adiabatic reversals 
(OSCAR) technique appears to be the shortest way to 
a single spin detection in magnetic resonance force mi- 
croscopy (MRFM) [l| . In this technique the cantilever tip 
(CT) vibrations in combination with a radio-frequency 
(rfj resonance field causes adiabatic reversals of the ef- 
fective magnetic field. Spins of the sample follow the 
effective magnetic field causing a small shift of the CT 
frequency, which is to be measured. The quasiclassical 
theory of the OSCAR technique has been developed in 
• The quantum theory of a single-spin measurement 
in OSCAR MRFM was presented in 7. 8]. 

It was shown in that the CT frequency can take 
two values corresponding to two possible directions of the 
spin relative to the direction of the effective magnetic 
field, B e f f . If the spin points initially in (or opposite to) 
the direction of B e ff, then the CT has a definite trajec- 
tory with the corresponding positive (negative) frequency 
shift. In the general case, the Schrodinger dynamics de- 
scribes a superposition of two possible trajectories - the 
Schrodinger cat state. It was also shown in that the 
interaction between the CT and the environment causes 
two effects. The first one is a rapid decoherence: the 
Schorodinger cat state transforms into a statistical mix- 
ture of two possible trajectories. Physically this means 
that the CT quickly selects one of two possible trajec- 
tories, even if initially the spin does not have a definite 
direction relative to the direction of B e ff- The second 
effect is an ordinary thermal diffusion of the CT trajec- 
tory. 

What was not taken into account in 00 was the effect 
of the direct interaction between the spin and the envi- 
ronment. In general, for an MRFM system one should 
consider two environments: one for the CT and the other 
for the spin. If the initial spin wave function describes 
a superposition of two spin directions relative to B e ff, 
then the spin generates two trajectories for the CT. The 
cantilever is a quasiclassical device which measures the 
spin projection relative to the direction of B e ff. The 



CT environment collapses the CT-spin wave function se- 
lecting only one CT trajectory and a definite direction 
of the spin relative to B e f f . This situation is similar to 
the Stern-Gerlach effect, but for a time-dependent B e ff- 
The direct interaction of the spin with its environment is 
extremely weak in comparison to the interaction between 
the CT and its environment. However, this weak interac- 
tion causes a deviation of the spin from the direction of 
B e f f after a collapse of the wave function. In turn, this 
deviation generates two CT trajectories. This scenario 
occurs again and again in the OSCAR MRFM. Normally, 
a collapse "forces" the spin back to its initial (after the 
previous collapse) direction relative to B e //- But some- 
times this collapse pushes the spin to the opposite di- 
rection, revealing the quantum jump - a sharp change 
of the spin direction and CT trajectory. It was shown 
in that the main source of the magnetic noise on 

the spin was associated with the cantilever modes whose 
frequencies were close to the Rabi frequency. 

There are two basic problems associated with the single 
spin OSCAR MRFM. The first problem is the theoreti- 
cal description of the statistical properties of quantum 
jumps. Unfortunately, the direct simulation of quan- 
tum jumps consumes too much computer time to be 
implemented. Recently, we have considered a simpli- 
fied model which describes the statistical properties of 
quantum jumps 9]. The second problem is more sophis- 
ticated: What is the characteristic time interval between 
two consecutive collapses of the the wave function? This 
Letter discusses the second problem. 

The CT position and momentum have finite quantum 
uncertainties. Thus, when the spin direction deviates 
from the direction of B e ff, the collapse does not occur 
instantly. During a finite time interval, the spin and the 
CT are entangled. The spin does not have a definite 
direction, and the CT does not have a definite trajec- 
tory. The dynamics of the CT-spin system on the time 
scale less than the time interval between two consecu- 
tive collapses can be described by the Schrodinger equa- 
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FIG. 1: Cantilever set-up. B ex t is the external static mag- 
netic field, Bi the rotating magnetic field, in the magnetic 
moment of the ferromagnetic particle glued on the cantilever 
tip, S the spin to be measured. 



tion. During this time, the average CT frequency shift 
is expected to be smaller (in absolute value) than the 
frequency shift corresponding to the definite direction of 
the spin relative to B e ff. We show that the experimental 
study of this effect can reveal the answer to a fundamen- 
tal problem of quantum dynamics: At what time does the 
collapse of the wave function occur if the quasiclassical 
trajectories are not well separated? 

Indeed, if the quasiclassical trajectories are initially 
well separated (the Schrodinger cat state) then the char- 
acteristic time of the collapse is the decoherence time, 
i.e. the time of vanishing of the non-diagonal peaks 
of the density matrix (the non-diagonal peaks describe 
the quantum correlation between two trajectories when 
these trajectories coexist during the same time interval 
fTotlTTj t. However, the Schrodinger cat state is a specific 
bizarre phenomenon in the macroscopic world. Namely, 
in a typical situation the collapse of the wave function 
occurs long before a well defined separation develops be- 
tween the two quasiclassical trajectories. There are a few 
simple cases, for which the exact solutions of the mas- 
ter equation have been obtained (see, for example, 
The exact solution describes a generation of the two qua- 
siclassical trajectories, their decoherence, and a thermal 
diffusion. Before the two quasiclassical trajectories are 
well separated, the exact solution describes the compli- 
cated dynamics of the density matrix elements. It is not 
clear if the master equation is capable to describe the 
wave function collapse when the two trajectories are not 
well separated. Even if the collapse time for this case is 
"hidden" in the solution of the master equation, we still 
do not know how to extract it analytically and numeri- 
cally. Only experiments could resolve this fundamental 
problem. We show that OSCAR MRFM could become 
one of these experiments. 

The OSCAR MRFM setup considered here is shown in 
Fig. ID 

The dimensionless Hamiltonian for the OSCAR 
MRFM in the rotating system of coordinates is taken 



as 

H = (p 2 x +x 2 )/2 + eS x + 2nxS z , (1) 

where we use the quantum units of the momentum 
Pq = h/Xo and the coordinate Xq — ^/hu} c /k c - Here oj c 
is the cantilever frequency, k c is the cantilever effective 
spring constant, e = ^B\/u) c , 7 is the magnitude of the 
electron gyromagnetic ratio and rj = (l/2)( r yX /u) c )G, 
where G = dB z /dx, is the magnetic field gradient at the 
spin location. The first term in describes the unper- 
turbed oscillations of the CT, the second term describes 
the interaction between the spin and the resonant rf field, 
and the last term describes the interaction between the 
spin and the CT. To take into consideration the magnetic 
noise on the spin caused by the thermal CT vibrations 
we add a random term A(t)S z to the Hamiltonian JTJ. 

The wave function of the system is assumed to have 
the form 

* = u a (x, t)ol + up(x, r)/3, (2) 

where a and (3 are the basis spin functions in the S z 
representation corresponding to the values S z = ±1/2, 
and t = uj c t is the dimensionless time. The Schrodinger 
equation splits into two coupled equations for u a (x, r) 
and up(x,T). If these two functions are identical (up to 
a constant factor) then the wave function can be repre- 
sented by a product of the CT and spin wave functions. 
In this case the average spin (S) has a magnitude 1/2. In 
the general case, the spin is entangled with the CT, and 
the average spin is smaller than 1/2. In our estimates we 
will use the following parameters from experiment pj: 

w c /2tt = 6.6kHz, k c = 600fiN/m, B 1 = 300^T, 

G = 4.3 x 10 5 T/m, X m = lOnm, T = 200mK, (3) 

where X m is the amplitude of the CT. Using these val- 
ues, we obtain the following values of parameters for our 
model 

X = 85/m, P = 1.2 x l(T 21 iVs, 77 = 0.078, , , 
x m = X m /X = 1.2 x 10 5 . (4j 

The relative frequency shift 5u>o = |Aw c /w c |of the can- 
tilever vibrations can be estimated to be 0, : 

Suj = 2Si] 2 /{2n 2 x 2 m + e 2 ) 1/2 = 4.7 x 10~ 7 . (5) 

(We insert the factor 2S for future discussion.) 

Suppose that initially the spin points opposite to the 
direction of B e f t . The frequency shift of the CT vibra- 
tions is then — Sluq. The magnetic noise acting on the spin 
causes a deviation of the spin direction from the direction 
of the effective magnetic field. Thus, it produces two tra- 
jectories of the CT with the frequencies ±8(Jq. Because 
of the quantum uncertainty of the CT position during 
the finite time (which we call the collapse time t c0 «), the 
wave function of the CT describes a single peak with an 



absolute value of the frequency shift less than Scuq . The 
collapse of the wave function changes the frequency shift 
to the value — 6u>o (or, sometimes, Sloq, in the case of a 
quantum jump). 

We simulate the quantum dynamics between two con- 
secutive collapses of the wave function of the system. In 
our simplified model the magnetic noise A(t) takes con- 
secutively two values, ±Ao. The time interval between 
two consecutive "kicks" of the function A(i) was taken 
randomly from the interval (— 3tr/4, 5tr/4), where tr = 
2ir/e is the dimensionless Rabi period. (In a more ad- 
vanced theory the characteristics of the magnetic noise 
should be derived from the parameters of the thermal CT 
vibrations.) The initial wave function is assumed to be a 
direct product of the CT and spin wave functions. The 
initial state of the CT is a coherent state 

l«o) =7r 1/4 ^2"/ 2 a "i7„( 2 :)exp{-(.T 2 + |ao| 2 )/2}, (6) 

n 

where H n (x) is a Hermite polynomial, and a$ = -4- (xq + 
Po); here Xq andpo are the quantum mechanical averages 
of x and p at r = 0. The initial direction of the spin 
is taken to be opposite to the direction of the effective 
magnetic field -B e // = is + 2h^xo, where i and k are the 
unit vectors in the positive x- and z-directions. 

In our numerical simulations we expand the functions 
u a (x,r) and u@{x,t) in (J2J over 400 eigenfunctions of 
the unperturbed oscillator Hamiltonian. During the time 
interval between two consecutive "kicks" of the noise 
function A(f) we have a time-independent Hamiltonian. 
Thus, we find the evolution of the wave function by diag- 
onalizing the 800 x 800 matrix and taking into considera- 
tion the initial conditions after each "kick" . The output 
of our simulations is the time interval At 3 = Tj — Tj—i 
between two consecutive returns to the origin for the av- 
erage value (x): {x{rj-\)) = and (x(rj)) = 0. To save 
computational time, we have used the values of param- 
eters e = 10, r] = 0.3, po — 0, x — 13. These values 
provide a relative large frequency shift of the CT oscilla- 
tions 0,@ 

6lu ps 7.9 x 10~ 3 . (7) 

The results of our simulations are shown in Fig. [21 
which demonstrates the deviation of Atj from the un- 
perturbed half-period of the CT oscillations it. We intro- 
duce Stj = \Atj — 7r|. With no magnetic noise (A = 0) 
the deviation 6Tj does not change with time 

Stj = St q = tt/Sluo ps 0.025. (8) 

If Aq 0, the value of 6tj decreases with time until the 
collapse of the wave function destroys the Schrodinger cat 
state. (Fig. 2 demonstrates the case for which the time 
interval between the collapses equals six half-periods of 
the CT vibrations). 

We will define the "effective frequency" u>j for each 
half-period of the CT vibrations: uij — tt/Atj. Next, 
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FIG. 2: Deviation of At, from the unperturbed half period 
of the CT oscillations it, as a function of the number of half 
periods j, for different Ao as indicated in the legend. Solid 
lines are the standard linear fits. On the vertical axis Srj = 
\A Tj -ir\. 

we introduce the effective frequency shift Sujj = tv/Stj. 
From Fig. 2 we can find the average effective frequency 
shift (Su>) between two collapses of the wave function. 

If our model of the noise were the adequate one, then 
from the experimentally measured quantity (Su>) we could 
determine the time interval between two collapses t co u. 
(It is clear from Fig. 2 that (Sui) is directly related to the 
value of T co ii for the assumed magnetic noise parameters.) 
Certainly, in a real situation this opportunity does exist 
if the average frequency shift (6u>) is significantly smaller 
than the expected shift Suq. We will call such situation 
"the case of the strong noise" . 

Note, that a decrease of the frequency shift may be 
interpreted as an effective decrease of the spin SS. Using 
Eq.|J5} we obtain 

6S = «<M - 5u; )(2^x 2 m + e 2 ) l ' 2 /{2 v 2 ). (9) 

Previously we have shown that in the case of strong 
noise, an experimentalist could determine the time in- 
terval T C oii between two consecutive collapses of the wave 
function by measuring the decrease of the frequency shift 
of the CT vibrations. Here we propose a special experi- 
ment which allows one to determine the time t co u for the 
case of weak noise in which the average frequency shift 
between two collapses is close to Suo- 

The problem is the following. Even a very weak noise 
generates a second CT trajectory with a frequency shift 
opposite to that of the first trajectory. Thus, two trajec- 
tories tend to separate at the same rate for any magnetic 
noise. Correspondingly, the time interval between two 
collapses is expected to be approximately the same for 
any level of magnetic noise. However, in the case of weak 
noise, the probability of the second trajectory is small, so 
its contribution to the average frequency shift becomes 
negligible. 

To overcome this obstacle, we propose an artificial 
change of the frequency shift using the "interrupted OS- 
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CAR technique" recently implemented in 0. In p|, the 
rf field is turned off for a time interval equal to half of the 
CT period, which is equivalent to the application of the 
7r-pulse which changes the direction of the spin relative 
to the effective magnetic field. We propose to turn off the 
rf field for the duration of the quarter of the CT period, 
which is equivalent to the application of the 7r/2-pulse. 
Suppose that initially the spin is parallel to the effective 
magnetic field, B e ff. If we apply a "7r/2-pulse" , the spin 
will become perpendicular to B e ff. Thus, we have two 
CT trajectories, each with the same probability. Before 
these two trajectories are separated, the CT will oscil- 
late with the unperturbed frequency, which is equal to 
one in our dimensionless units. After the collapse the 
frequency shift is ±Su>o, with equal probabilities. Thus, 
using a "7r/2-pulse" we can achieve a maximum possible 
reduction of the frequency shift. If we apply a periodic 
sequence of "7r/2-pulses" with the period r p (r p > T co u), 
then the average frequency shift (Slo) is 

(6u>) = 5w (t p - t coU )/t p . (10) 

Manipulating t p one can achieve a significant decrease 
of (Sui) in comparison with Suq. Using Eq. fTT)|) one 
can determine the collapse time from the experimental 
value (Slo). Thus, the collapse time can be measured 
for the case of weak magnetic noise. Finally, based on 
quasiclassical theory [j| we will estimate the reduction 
of the average frequency shift caused by the noise for 
the experimental conditions 0. The amplitude of the 
thermal CT vibrations near the Rabi frequency can be 
estimated to be: 

a T = {u c /u a ){k B T/2k ) 1 ^ = 38/m. (11) 

The square of the characteristic spin deviation during a 
single reversal (half of the CT period) is 

(A6»i) 2 ps ZAGa 2 T /{uj c X m ) = 9 x 1CT 7 . (12) 

We have no idea about the order of the collapse time t co u . 
If we assume that the collapse occurs when the separation 
between the two trajectories with the frequencies 1 ± 
Sojq is equal to 1/2 (the quantum uncertainty of the CT 
position in the coherent state is ((Ax) 2 } = 1/2), then we 
obtain for t co u 

t coU shiT coU ~ 1/ (4x m Su> ) ps 4.4. (13) 

It follows from Eq. (|13fl that the wave function collapses 
during the second period of the CT vibrations. If we 
estimate the probabilities of the two CT trajectories as 



Pi ~ 1 — (A6>i) 2 (for the trajectory with the initial fre- 
quency shift) and P2 ~ (A6*i) 2 (for the trajectory with 
the opposite frequency shift), then the average CT fre- 
quency shift can be estimated to be 

(Sw) = <MPi - P2) = *w (l - 2(A^) 2 ). (14) 

This estimated reduction of the CT frequency shift is 
clearly negligible. 

Consider the opposite extreme case. Suppose that the 
collapse occurs when the separation between the two tra- 
jectories is of the order of the thermal CT fluctuations 
(ksT/kc) 1 ^ 2 ps 150pm or 1760 in dimensionless units (as 
before, we used the values of parameters in (3)). In this 
case, the time interval between the two consecutive col- 
lapses t co u is of the order of 10 4 periods of the CT oscil- 
lations. During this time, the characteristic spin devia- 
tion is (AO) 2 ~ 0.02. Then, we have for the probability 
Pi ps (A#) 2 /4 ps 5 x 10~ 3 . The average frequency shift 
is (6u>) ps 6ljq(1 — 10~ 2 ). One can see that even in this 
extreme case the reduction of the frequency shift is ex- 
pected to be small. 

Thus, the experimental conditions in p] probably cor- 
respond to the case of the weak noise. In such a situation 
the collapse time could be measured using the periodic 
sequence of "7r/2-pulses" described earlier. 

We have demonstrated a procedure for measuring the 
mysterious collapse time in the OSCAR MRFM tech- 
nique. We simulated the quantum dynamics of the spin- 
CT system. Unlike the previous studies of the quantum 
dynamics we took into consideration the direct interac- 
tion between the spin and the environment (the mag- 
netic noise). This noise causes (i) a deviation of the 
spin from the direction of the effective magnetic field 
and (ii) entanglement between the spin and the CT. For 
the case of weak magnetic noise, the same effect can be 
achieved the effective "7r/2-pulses" with the "interrupted 
OSCAR" technique. The spin-CT entanglement influ- 
ences the frequency of the CT oscillations before the wave 
function collapse takes place. This effect can be described 
as an effective decrease of the single spin magnitude. We 
demonstrated that the experimental measurement of the 
OSCAR MRFM frequency shift could reveal information 
about the time interval between two consecutive collapses 
of the wave function. 
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